Single to Double Hump Transition in the Equilibrium Distribution Function of 

Relativistic Particles 
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We unveil a transition from single peaked to bimodal velocity distribution in a relativistic fluid 
under increasing temperature, in contrast with a non-relativistic gas, where only a monotonic broad- 
ening of the bell-shaped distribution is observed. Such transition results from the interplay between 
the raise in thermal energy and the constraint of maximum velocity imposed by the speed of light. 
We study the Bose-Einstein, the Fermi-Dirac, and the Maxwell- Jiittner distributions, all exhibiting 
the same qualitative behavior. We characterize the nature of the transition in the framework of 
critical phenomena and show that it is either continuous or discontinuous, depending on the group 
velocity. We analyze the transition in one, two, and three dimensions, with special emphasis on 
two-dimensions, for which a possible experiment in graphene, based on the measurement of the 
Johnson-Nyquist noise, is proposed. 
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Back in 1911, F. Jiittner derived a relativistic analogue 
of the Maxwell-Boltzmann equilibrium distribution for 
classical (non-relativistic) gases. To this purpose, he re- 
sorted to an entropy minimization procedure, subject to 
the relativistic energy-momentum constraints [l|- This 
top-down (macro-to- micro) derivation, left room for some 
debate on whether the exact form of the equilibrium dis- 
tribution for relativistic particles was the one proposed 
by Jiittner, or rather some variant thereof 0, Very 
recently, conclusive evidence for the original form pro- 
posed by Jiittner has been brought by numerical sim- 
ulations of fully relativistic molecular dynamics in one 
and two dimensions 0, [H[ . This is all but an academic 
exercise, since the Jiittner distribution is known to play 
a major role in the interpretation of current and future 
experiments in many sectors of modern physics, such as 
quark-gluon plasmas produced in heavy- ion collisions @ , 
relativistic astrophysics Q , distortions of the cosmic mi- 
crowave background [U, and lately, possibly also in the 
study of electron flows in graphene @, [Io| . 

In this Letter, we wish to call the attention on an ap- 
parently hitherto unexplored property of the Jiittner dis- 
tribution, namely the fact that, under the constraint of 
an increasing temperature (ultra-relativistic limit mc 2 < 
ksT), the one-dimensional distribution develops a transi- 
tion in velocity space, with the emergence of two separate 
peaks, moving with opposite, non zero-speeds. For two 
and three dimensions, it generates, in velocity space, a 
ring and a hollow sphere (see Fig. [TJ), respectively. This 
stands in sharp contrast with the way a non-relativistic 
gas at rest responds to the constraint of an increasing 
temperature, namely through a progressive broadening 
of the Gaussian shape, which enhances the high-speed 




FIG. 1. Three-dimensional Maxwell- Jiittner velocity distri- 
butions, according to Eq. JTJ, for two different temperatures. 
The blue (£ = 20) and gray (£ = 1) isosurfaces stand at 1/3 
and 1/4 of their respective maxima. Shown in the inner region 
of each isosurface, is the color gradient of the distribution, 
with red and blue colors denoting high and low concentration 
of particles, respectively. Note that the spheres are plotted in 
velocity space, so that the maximum radius for the external 
isosurface is close to \v\ = 1. 



population, leaving nonetheless the least-energetic, zero- 
speed, particles as the most probable population, since 
this is the one best conforming to the zero net-motion 
constraint. On the contrary, the transition exhibited 
by the Jiittner distribution, signals that, even in a gas 
at rest, the temperature constraint can only be met by 
clustering most of the particles around two oppositely 
moving beams, thereby depleting the zero-speed parti- 
cles in the process. The above considerations readily 
generalize to the case of a moving gas, the main change 
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being that the two oppositely moving beams get differ- 
ently populated, the co-moving one being enhanced and 
the counter-moving being correspondingly depopulated. 
This phenomenon is quite general, and it might apply to 
a whole class of systems where physical signals are forced 
to move close to the their ultimate limiting speed. In the 
following, we provide mathematical details of this transi- 
tion and also discuss conditions under which it could be 
experimentally probed in graphene experiments. 

The probability distribution function of particle veloc- 
ities in d-dimcnsional rclativistic gas is described by the 
following single-particle distribution function (in natural 
units (c = H = e = k B = m = 1) [U E3: 
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where 7(1;) = l/vl — v 2 is the Lorentz factor, T the 
temperature, U the group velocity, v the velocity of the 
particles, A a normalization constant, v = \v\, and U = 
\U\. The subscript A denotes the Fermi-Dirac (A = 1), 
the Bose-Einstein (A = —1), and the Maxwell- Jiittner 
(A = 0) distributions, respectively. At the moment, we 
will neglect the chemical potential (fx = 0). The ther- 
mal behavior of a rclativistic gas at equilibrium is best 
characterized by the parameter £ = mc 2 /ksT = l/T, 
that is commonly used to differentiate between the ultra- 
relativistic (£ S> 1) and the relativistic (£ <C 1) regimes. 

Under the constraint of a limiting velocity imposed by 
the theory of relativity, the three distributions share an 
interesting property in the temperature dependence. In 
particular, above a critical temperature (T c = l/£ c ), the 
shape of the function changes from a nearly Gaussian to 
bimodal (see Fig. [5]). As a result, while below the critical 
temperature the majority of the particles move at speeds 
close to the group velocity U (zero in Fig. [2]) , above crit- 
icality two populations of particles emerge, with a veloc- 
ity distribution sharply peaked around opposite speeds, 
close to the speed of light. Here, we show that this change 
in the functional dependence of the velocity distribution 
can be described in the framework of a transition, which 
might be either discontinuous (first order) or continuous 
(higher order), depending on the group velocity U. In 
two and three dimensions, the same qualitative behavior 
is observed, where, instead of two peaks, a ring (in two 
dimensions) and a hollow sphere (in three dimensions, 
sec Fig. [T]) is obtained at higher temperatures. In both 
cases, at sufficiently high temperatures, the radius in ve- 
locity space reaches the value corresponding to the speed 
of light (see Fig. [[]). 

To characterize the transition, we introduce an order 
parameter defined as the distance A between the peaks, 
such that A = in the single peaked distribution and 
A ^ in the bimodal one. Let us begin by consider- 
ing the case U = 0, and measure the order parameter 
dependence on the temperature, T = l/£, as shown in 
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FIG. 2. Bose-Einstein (A = -1), Maxwell- Jiittner (A = 0), 
and Fermi-Dirac (A = 1) velocity distributions, for d — 1, 
according to Eq. (fjl. The distributions correspond to a rel- 
ativistic fluid, with group velocity U = and two different 
temperatures (T = l/£), namely, £ = 1 and £ = 20. 



Fig. [31(a). From this figure, we can appreciate that, be- 
low the critical temperature T c = l/£ c , the distribution 
function has only one peak, A = 0, while above T c , a bi- 
modal profile develops, with the order parameter growing 
from zero to an asymptotic value A = 2. In this limit, 
the width of the distribution shrinks to zero and the dis- 
tribution is a superposition of two Dirac deltas, literally 
corresponding to a discrete fluid moving at ±c. 

In the inset of Fig. Efa), we analyze the singularity 
at T = T c , namely, we plot the order parameter depen- 
dence on the rescaled control parameter l/£ — l/£ c . A 
continuous transition is observed, with the order param- 
eter being zero at the critical temperature and growing 
according to a power law A ~ (l/£ — l/£ c ) ' 5 above it. 
This exponent corresponds to the inverse of the expo- 
nent 5 in the theory of critical phenomena [l3j . The 
same qualitative behavior is observed in two and three 
dimensions. Below, we describe the way the exponent 
and critical temperature can be obtained analytically, in 
the limit U 0. 

Due to the fact that the aforementioned transition is 
driven by the j d+2 pre-factor in the distribution, which 
is symmetric around v = in velocity space, one can, 
without loss of generality, calculate A (diameter of the 
ring, for d = 2, or of the hollow sphere, for d = 3) along 
the direction v x . In the limit U = 0, from the calculation 
of the maxima of the distribution, Eq.([T]), one obtains 
the trivial solution v x = and two additional ones, cor- 
responding to the solutions of the algebraic equation, 
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non-zero solutions, wc obtain 



y/[2 + d + W(q)]*T2-l 



(3) 



T(2 + d+W{q)) 

where q = (2 + d)e~^ 2+d 'A and W(x) is the Lambert 
W-function. From this equation, we can also obtain the 
critical temperature as, 

1 

Tc ~ d + 2 + W(q) ' (4) 

where, in the case of the Maxwell- Jiittner distribution 
(A = 0), i.e., q = W(q) = 0, and so T c = l/(d + 2). 
In general, for this distribution, the temperature depen- 
dence of A is given by, 

A = -^-^s(e + 2T c ) , (5) 

& ~T J- C 



with e = T - T c . For T > T c , A - y/8e /T c , where the 
exponent is 1/2, in excellent agreement with the numeri- 
cal data in the inset of Fig. E£a). This exponent depends 
neither on the spatial dimension, nor on A, i.e., it is the 
same for the three distributions. However, the critical 
temperature T c depends on these quantities, and can be 
obtained analytically for each value of d and A. 

Figure [3Jb) shows the numerical results for the tem- 
perature dependence of A at different group velocities 
U, for the Maxwell- Jiittner distribution. For any U ^ 0, 
a jump on A is observed at the onset of the transition, 
resembling thermal first-order transitions. In Fig. 21 we 
plot the size of the jump Aj as a function of U , where a 
power law is obtained Aj ~ , with (3 = 0.32 ± 0.02. 
This exponent corresponds to the order parameter expo- 
nent P in the theory of critical phenomena 
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the error bars, the same exponent was obtained for all 
combinations of A and d, suggesting that (3 is indepen- 
dent of these two parameters. Nevertheless, as shown 
in the inset of Fig. EJb), the transition temperature in- 
creases with U and decreases with d. 

The above transition is not transmitted to the con- 
served macroscopic quantities (energy-momentum) and 
consequently, it is not straightforwardly observed at the 
macroscopic level. Instead, a microscopic analysis is re- 
quired. In order to experimentally detect the transition, 
thermal energies of the order of the rest energy of the 
gas particles (£ ~ 1) would be required. This could 
be achieved, for example, in hot electron plasmas. In 
particular, it is of great interest to explore whether this 
transition can contribute to a deeper understanding of 
small condensed matter systems with potential techno- 
logical applications. A good candidate in this respect 
is graphene. Since its discovery [l4|, 15 1, graphene has 
continued to surprise scientists with an amazing series 
of spectacular properties, such as ultra-high electrical 
conductivity, ultra-low viscosity to entropy ratio, com- 
bination of exceptional structural strength and mechan- 
ical flexibility, and optical transparency. It consists of 
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FIG. 3. (Top panel) Temperature (T = l/£) dependence of 
the order parameter A, defined as the distance between peaks, 
for the Maxwell- Jiittner distribution (A = 0), in different spa- 
tial dimensions d and group velocity U = 0. Shown in the 
inset is a double-logarithmic plot of A as a function of the 
rescaled control parameter l/£ — l/£ c , where £o stands for 
the transition temperature. (Bottom panel) Temperature de- 
pendence of the A for different group velocities. The inset 
shows the transition temperature T c = l/£ c dependence on 
the group velocity U. 



literally a single carbon monolayer and represents the 
first instance of a truly two-dimensional material (the 
"ultimate flatland" [HI ) , where electrons move like mass- 
less chiral particles and their dynamics is governed by 
the Dirac equation, following the dispersion relation, 
E{k) = shvp \k\, where k = (k x ,k y ) is the wave vector. 
The constant s = ±1 distinguishes between electrons (+) 
and holes (— ), and vf is the Fermi velocity that plays the 
same role in graphene as the speed of light in relativity. 
This relation implies that carriers always move at the 
same Fermi speed, regardless of the Fermi energy. For 
simplicity, we will work only with the electronic density 
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FIG. 4. Size of the jump Aj as a function of the group velocity 
U for different dimensions d. The quantity denotes the slope 
of the best fit. Shown in the inset, is the two-dimensional 
Fermi-Dirac distribution according to Eq. Q for £ = 61.6 
(/i-BN graphene at 5K). 



(s = 1), the extension to include holes being straightfor- 
ward. 

In our context, pristine graphene corresponds to the 
ultra-relativistic limit, where the velocity distribution 
function consists of two Dirac deltas at ±i>f, so that 
the transition cannot be observed. However, the elec- 
tronic spectrum of graphene changes depending on the 
substrate. For example, on SiC the energy spectrum 
presents a gap of width 2mv F = 0.26cV and on /i-BN 
(hexagonal boron nitride) a gap of 53meV (l7l - fl9j . and 
can be manipulated by constructing graphene nanorib- 
bons, where the energy gap depends on the width of the 
ribbon [UELI- 

A gap in the spectrum of graphene corresponds to non- 
zero mass of the electrons, so that the dispersion rela- 

a 2 , where a = wivf/H. 



tion becomes E a {k) = hvp\[W 
With this dispersion relation, the velocity of the electrons 
can be calculated as v = (l/h)dE(k)/dk, thus obtaining, 

v = vpk/\/ k 2 + a 2 . Note that the velocity is no longer 
constant and depends on the Fermi energy via k. The 
electronic density n is defined by 



(6) 



where f(k) = [1 + exp((E(k) - ^)/fc s T)]- 1 is the Fermi- 
Dirac distribution. To calculate the density of states in 
velocity space, we change variables, from the wave to 
the velocity vector spaces, obtaining n = J D(v)f(v)d 2 v, 
where D(v) = a 2 j 4 / '(2tt) 2 vf is defined as the density 
of states (DoS). The Fermi-Dirac distribution function 



takes the form, 
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with rj = n/kBT, and the parameter £ now defined as 
£ = mVp/ksT. Note that the DoS tends to push all the 
particles moving at the Fermi velocity, while the Fermi- 
Dirac distribution counters this effect, given the limit 
imposed by the Fermi speed. At low temperature, the 
chemic al p otential can be approximated by the Fermi 
energy [22j . Ep, which can be tuned experimentally by a 
gate voltage for low concentration of electrons [23| . Since 
low temperatures are required, to avoid electron-phonon 
interactions, we characterize the transition of the velocity 
distribution functions by changing the Fermi energy, i.e. 
7], instead of the temperature. With this change in the 
control parameter we also have the advantage that the 
speed at which the Fermi-Dirac distribution attains its 
maximum corresponds, to a good approximation, to the 
maximum speed at which carriers can move (see inset of 
Fig. QJ. Therefore, we expect that the transition can be 
observed by measuring the thermal or Johnson-Nyquist 
noise [24H26| . Let us assume that we have a typical sam- 
ple of graphene on /i-BN, and the electronic density is 
manipulated with an external gate voltage by the rela- 
tion n ~ kt] [27]], where k is a proportionality constant 
that depends on the electric capacitance of the substrate 
and the temperature. Thus, for low carrier concentra- 
tion, even if there is no drain voltage, there are current 
fluctuations SI around zero due to the thermal motion 
of the electrons in the sample. The maximum amplitude 



for these fluctuations can be written as SI„ 



env„ 



where e and I are the electric charge of the electrons and 
the cross section of the sample, respectively. We define 
the dimensionless maximum amplitude for the current 
fluctuations as T = 25I max / enlvF- 

For a fixed low temperature, we can observe from Fig. [5] 
that there is a critical rj c above which T emerges. This 
critical value increases with decreasing temperature. We 
have considered different temperatures, T = 20,30,50 
K. Note that the critical rj c = 1.05£ + 0.26, is almost 
the same as £, by a proportionality constant 1.05, and 
the residual is just due to the non-linear behavior of the 
curve, close to the critical temperature where, due to the 
thermal energy, the thermal noise appears, regardless of 
77. Note that this expression is tantamount to stating 
that the Fermi level must be higher than half of the gap 
energy in order to have electrons in the conduction band. 
The exponent of the continuous transition is also 0.5, and 
is independent of the temperature of the sample in the 
regime of low temperatures. 

Some transport properties in graphene may also de- 
pend on the velocity instead of the momentum, e.g. sat- 
uration current in ballistic regime 28|, 29| and mobility 
of carriers [3(|. Therefore, it is plausible to expect that 
the transition in the Fermi-Dirac distribution discussed 
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FIG. 5. Dimensionless maximum amplitude for the current 
fluctuations T for gaped graphene on a substrate A-BN as 
a function of n. The inset shows the critical n c , for which 
the thermal noise emerges, as a function of the gap-induced 
mass, described by £. The solid line represents the equation 
n c = 1.05| + 0.26. 



in this Letter, might play an important role in developing 
a better understanding of (some of) the spectacular prop- 
erties of graphene. Furthermore, with the aim of build- 
ing graphene transistors at nano-size scales, the study 
of nanoribbons (GNR) [2^, 21 1 has become very popu- 



lar and, therefore, the individual carrier dynamics might 
play a major role in affecting the transport properties of 
these devices. 

Summarizing, we have shown that the equilibrium dis- 
tribution for relativistic particles, no matter whether 
classical or quantum, exhibits a transition as the temper- 
ature is brought close to the rest energy (ultra-relativistic 
limit). This transition is the organized response of the 
distribution of particles to the constraint of an increas- 
ing temperature, compatibly with the existence of a lim- 
iting speed for the propagation of physical signals. We 
have also discussed conditions under which such transi- 
tion could potentially be detected in current and future 
graphene experiments. In the case of graphene, we have 
found that the transition takes place, not only by increas- 
ing the temperature, but also by increasing the Fermi 
energy at costant low temperature. 
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